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We discuss how the simultaneous crossovers of deconfinement and chiral restoration can be real- 
ized. We propose a dynamical mechanism assuming that the effective potential gives a finite value 
of the chiral condensate if the Polyakov loop vanishes. Using a simple model, we demonstrate that 
our idea works well for small quark mass, though there should be further constraints to reach the 
perfect locking of two phenomena. 
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' Introduction Chiral symmetry plays an important role in effective model approaches to Quantum Chromody- 
namics (QCD). The hadronic properties at low energy have been successfully described by chiral effective models 
fN| ■ such as the linear sigma modelW , the Nambu-Jona-Lasinio (N JL) model 0, , the chiral random matrix model 
chiral perturbation theory and so on. The nature of the chiral phase transition at finite temperature can be 
Q-c classified by chiral symmetry according to the universality argument and investigated in these effective models • 
1 In particular, in the massless two-flavor case, we can expect a chiral phase transition of second-order that belongs to 
the same universality class as the 3d 0(4) spin model. Then, we can anticipate what could occur in the real world 
with finite but sufficiently small (up and down) quark masses. 

The deconfinement phase transition is rather obscure and veiled because it is only well-defined in the heavy quark 
■ limit, which is too far from the real world. In the heavy quark limit, i.e., in the absence of dynamical quarks, 
£> \ the Polyakov loop serves as an order parameter for deconfinement and the phase transition is characterized by the 
spontaneous breaking of center symmetry . In the presence of dynamical quarks the center symmetry is explicitly 
broken and no order parameter or criterion has been established for the deconfinement transition As mentioned 
above, on the other hand, the chiral phase transition at high temperature or baryon density has been well-understood 
by means of effective models. Those model studies based on chiral symmetry, however, lack any dynamics coming 
, from the Polyakov loop , except for some efforts to clarify the interplay between chiral dynamics and the Polyakov 

loop pi mill mm. 

It is often argued that the mixing between the Polyakov loop, L, and the chiral order parameter (chiral condensate), 
X, can account for the lattice QCD observation that deconfinement and chiral restoration occur at the same pseudo- 
critical temperature |l5j| . The mixing argument is, however, not sufficient to give a satisfactory explanation on the 
Oh, lattice QCD data. In the (T, m q ) plane (m q being the current quark mass), as discussed in jl6j| . there appear two 
terminal points of the first-order phase boundary, namely, the critical end-points (CEPs). One is the chiral CEP 
denoted by (Tq, m q ) and the other is the deconfinement CEP denoted by (Tq, nig). For example, in the two-flavor 
three-color case, it is known that (T c - 170 MeV, m c q = 0) and (T D - 270 MeV, m q ~ I GeV) 0[l7|]. It should be 
noted that the CEP is a true second-order critical point with a divergent susceptibility. Mixing means that L and x 
should share the same singularity in their susceptibilities near the CEP; for m q = m q (and m® ), the susceptibilities 
of L and x both diverge at T — Tq (and To respectively). The susceptibility peak is smeared as m q leaves from the 
CEP. The important point is that, with m q fixed at a certain value near the CEP, a smeared bump originating from 
the other CEP may be observed separately as well as a sharp peak from the closer CEP. Thus the mixing argument 
cannot exclude a double-peak structure. For m q ~ m q for example, the Polyakov loop susceptibility can have a sharp 
peak around Tq (coming from the mixing with the diverging chiral susceptibility) as well as a broad bump around 
Tq (coming from a remnant of deconfinement). We should be cautious about the mixing argument to understand the 
lattice QCD data in which no double-peak structure has been seen for any m q [TEl lfl| . 

In fact, the locking between the two crossover phenomena depends on the detailed properties of the interaction. 
The purpose of this letter is to propose a simple mechanism to exclude the undesirable double-peak structure and to 
complement the shortcomings of the mixing argument. 

Idea To make our idea clear in general setting, let us suppose that we have a full effective potential, i.e., 
V c s[L,x] m q]- I n principle, the behavior of L and x should be completely determined by V c g [L, x; m q ]. Thus the 
question arises; what property of 14ff[L, x; rn q ] can give rise to the simultaneous crossovers? Our idea is as follows. 

First of all, an important property follows from the theoretical arguments given by Casher and 't Hooft [l9| . 
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According to their arguments, the confined phase must have a non- vanishing chiral condensate, which suggests that the 
chiral phase transition should occur at higher temperature than deconfinement. This means that V e g [L = 0, x; m q = 0] 
leads to x 7^ at any temperature if L = is imposed by hand (or approximately chosen as a minimum of the effective 
potential). This property has not been proven in QCD (see also [20j) but is realized in the strong coupling analysis 
[13j and assumed here. 

Next, because L has turned out to behave approximately as an order parameter in lattice simulations [lfij . we 
can expect that L is almost zero below the deconfinement crossover temperature, Td, regardless of dynamical quarks. 
Then, together with the above property, x must have a non- vanishing value below Td even for m q ~ (i.e., spontaneous 
chiral symmetry breaking). Thus the chiral restoration temperature, T x , is greater than or equal to Td- 

In contrast, the critical temperatures at the CEPs are T D ~ 270 MeV (for m q = oo) > Tq ~ 170 MeV (for m q = 0), 
which implies Td > T x for an intermediate value of m q , unless chiral or center symmetry is overwhelmingly broken. 

Our idea is that Td = T x is likely to be realized by T x > Td from the properties of the effective potential and 
Td > T x from, so to speak, the boundary condition. Chiral symmetry is broken by m q ^ 0, while the center symmetry 
breaking is suppressed by the constituent quark mass even for small m q . Hence, our idea is expected to work especially 
for m q ~ 0. This mechanism can complement the mixing argument and lead to a robust single-peak structure in the 
susceptibilities. 

Effective Model For the purpose of demonstrating our idea, we propose a simple chiral effective model with 
Polyakov loop dynamics. If the Polyakov gauge (A4 is static and diagonal) is employed, the Polyakov loop (or strictly 
speaking its phase) appears in the quark action as an imaginary quark chemical potential |2'2t |23| . Thus we can 
uniquely determine the coupling between the Polyakov loop and quark excitations. We shall take this advantage by 
adopting the NJL model which is given in terms of quark degrees of freedom. 

The conventional Lagrangian density of the NJL model is 

Q 

£ N JL = q^Yd^ - m q )q+-{(qq) 2 + {qh 5 Tq) 2 }, (1) 

where m q = 5.5 MeV and G — 2 x 5.496 GeV~ 2 . The momentum integration is regulated by the cut-off A = 631 MeV. 
These model parameters are chosen as to reproduce the pion mass and decay constant at zero temperature Q . In the 
mean field approximation the thermodynamic potential is given by 

n N j L /V=^(M- mq f-2N c N f J ^{T p + Tln[l + e-^)/ T ]+Tln[l + e-^)/ T ]}. (2) 

V is the spatial volume, N{ is the flavor number fixed as Nt = 2 throughout this letter, and fi is the quark chemical 
potential. We neglect any \x dependence in G as usual Q . The energy of quasi-quarks is given by E p = VP 2 + A/ 2 
with the constituent quark mass M = m q — G(qq). The cut-off is imposed only on the first term in the curly brackets 
(zero-point energy) in the present analysis. The finite temperature contribution has a natural cut-off in itself specified 
by the temperature. Identifying the imaginary quark chemical potential with the Polyakov loop, we can define our 
model by the following thermodynamic potential, 



n/V = V gluc [L] + ^-(M- 



2G 

3„ 



\2 



2N r N f 



J ^3{T p + T^Tr c ln[l + Te-^)/ T ] + T^Tr c ln[l + it e -(B,+/0/T] ^ (3) 



{2-Kf 

where the Polyakov loop is an SU(A C ) matrix in color space explicitly given by 



L(x) — Texp 



-i / dx4 ^4.4(2:4, x) 
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This coupling between the Polyakov loop and the chiral condensate can be derived also in the strong coupling approach 
[Til ll3T|. The generalization to aQCD (QCD with dynamical quarks in the adjoint representation in color space) is 
readily available if one replaces the Polyakov loop by the Polyakov loop in the adjoint representation (c.f. |13)- I n 
this letter we will focus only on the case in the fundamental representation. 

Since the four-quark coupling constant, G, contains the information on gluons, G should depend on L. We simply 
neglect this possible L dependence. Nevertheless, this L dependence makes no qualitative difference because G 
incorporates all gluons and would not be much affected by L only, where L is essentially the temporal component of 
gluons. This approximation is acceptable in the same level as neglecting the possible fj, dependence in G. 

It is worth noting that the NJL model with both quarks and the Polyakov loop is not incompatible with confinement 
at low temperature. To make this clear, let us assume that confinement corresponds to the condition, L = 0, though 
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this is not precise due to the center symmetry breaking. The Taylor expansion of the logarithmic terms generates 
a series in powers of L whose power corresponds to the quark excitation number. L, L 2 ,L 4 ,... terms have non- 
trivial triality and break the center symmetry, while L 3 , L(— 1), . . . terms have zero triality (color-singlet part) 
corresponding to baryons, mesons, etc. In a more elaborated approximation scheme used in 0, Il3| . the group 
integration of L with a vanishing mean-field, L = 0, singles out only the terms with zero triality. [Our notation is 
sloppy as long as no confusion arises. L is used for both the Polyakov loop matrix and its traced expectation value.] 
Roughly speaking, the limit of L = allows only excitations with zero triality such as baryons and mesons. In the 
present treatment the same argument holds in part (see ©). Although our model © is described in terms of colored 
quarks, the thermal excitation must consist of colorless composites as long as L is small. Once L gets larger at higher 
temperature, our model would reduce into the standard NJL model allowing thermal quark excitations. 

The remaining part, V^i uc [£], is the effective potential only in terms of the Polyakov loop. Here we shall adopt a 
simple choice for V^i uc [L] at the sacrifice of quantitative accuracy; we employ the leading order result of the strong 
coupling expansion that is simple and yet reasonable qualitatively as compared with the lattice results [24j. This 
potential has only one parameter a (the lattice spacing). For _/V c = 3 we can write it as 



V giuc [L] ■ a 3 /T = -2(d- 1) c _<7a/T |Tr c L| 2 - In [- |Tr c L| 4 + 8Re(Tr c L) 3 - 18|Tr c L 
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with the string tension a — (425 MeV) 2 . The first term comes from the kinetic part and the second term is just 
the logarithm of the Haar measure associated with the SU(3) group integration. V^i uo [L] leads to a first order phase 
transition with the critical coupling 2(d — l)e~ cra / Td = 0.5153. We can fix the deconfinement transition temperature 
as the empirical value T d = 270 MeV by choosing a -1 = 272 MeV. 

Apparently, the present model has two cut-offs, i.e., A and a -1 . This means that the model has two independent 
scales for chiral symmetry breaking and confinement, or for mesons and the Polyakov loop. Since QCD has only one 
scale, A and a -1 should be related to each other in principle. Here we simply fix them as model parameters since we 
are working in the effective model to abstract the essence for each dynamics [2l| . 

Numerical Results Before we refer numerical calculations in the N c = 3 case, let us introduce an ansatz 
to simplify the analysis. In the Polyakov gauge we can parametrize the SU(3) Polyakov loop matrix as L = 
diag(e i<#, ,e i *',e- i ^ +<# '')). The perturbative vacuum has 4> = <\>' — and we can choose the confining vacuum at 
4> = 27r/3, 4>' = [2^. Thus we shall fix 4>' = from the beginning for simplicity. Once this ansatz is accepted, we 
can rewrite the potential © only in terms of the traced Polyakov loop, i.e. I = (Tr c L)/A c = (1 + 2cos0)/3. After 
straightforward calculations we can reach the expression; 



Tr r In 
= ln 



Tr r In 



l + £t e -(Sp+/0/r 

1 + (3 I - l)e-^-">/ T + e -2(E p -^)M + ln [! + ( 3 l _ 1)c -(E P +K)/T + e ~2(E p+ »)/T 

+ lnfl + e-lZp-MT] inh + c -(£ p +m)/t 



(6) 
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FIG. 2: The order parameter and susceptibility around the chiral critical end-point with /ie = 321 MeV. 



When 1 = 1, the model is reduced into the standard two-flavor NJL model having the chiral phase transition at 
T x = 175 MeV in the chiral limit. If I is forced to be zero by hand, the temperature effect is so suppressed that 
spontaneous chiral symmetry breaking can sustain until T x ~ 520 MeV, that is much higher than Tq. Therefore the 
chiral phase transition cannot occur until the Polyakov loop jumps from nearly zero to a certain finite value. Our 
model (01 satisfies the essential assumption in our idea, T x > Td, for m q ~ 0. 

Figure ^ (left) shows the resulting behavior of order parameters as functions of temperature at fj, = 0. To see 
the simultaneous crossovers clearly the chiral condensate x — (qQ.) i s normalized by the value at T = (denoted by 
Xo)- The results from fl are represented by the solid (for x/xo) and dashed (for I) curves. The dotted curves are 
the results from I^njl and ^giue without any interaction between x and I for reference. It should be noted that the 
mixing interaction between x and I vanishes at zero temperature in this model. Consequently the normalization xo 
is identical for both the chiral condensates from f2 and f^NJL- 

In the presence of dynamical quarks, as seen from the figure, the Polyakov loop shows a crossover around the 
pseudo-critical temperature T c ~ 200 MeV. At the same time the chiral condensate is affected by the Polyakov loop 
such that it tends to be almost constant as long as T < T c . The pseudo-critical temperature can be read from the peak 
position of each susceptibility. Here we shall define the dimensionless susceptibility of the chiral order parameter and 
the Polyakov loop by using the curvature inferred from the potential @. First, we define the dimensionless curvature 
matrix C by 



_ A 2 d 2 PSI 
qq ~ dM 2 A 3 !/' 



C u = 



d 2 0SI 
dl 2 A 3 V' 



C q l — Cla — 



Iq 



Ad 2 

dMdlA^V' 



Roughly speaking, C q ^ (C u 1 ) corresponds to the chiral (Polyakov loop) fluctuation and C q 
of a quark and the Polyakov loop. Then, the susceptibility is given by the inverse of C; 



(7) 



is the interaction vertex 



Xq 



C u 



CqqCu 



(I 2 ' 



Xl 



c„ 



CqqCu 



C 



2 ■ 



(8) 



The physical meaning of the above equations is transparent if we notice that the fraction can be expanded as Xq — 
Cqq 1 + C~ q 1 CqiC^ 1 CiqC~ q 1 + • • ■ , that is the sum over mixing contributions. The mixing pattern is similar to the 
relation between the chiral susceptibility and the baryon number susceptibility H^,|2^|. As shown in Fig. ^ (right), 
our idea works pretty well so that the peak of the Polyakov loop susceptibility can be found just near the peak of the 
chiral susceptibility, though perfect coincidence is not reached. 

When the quark chemical potential, [i, becomes larger, the Polyakov loop shows a crossover with smoother slope 
because the density effect itself breaks the center symmetry explicitly. It is widely accepted that the chiral phase 
transition becomes of first-order at large \i. In the (//, T) plane, therefore, we can expect another CEP |23, [23 • 
Actually we found the CEP at (/i£ = 321 MeV, Te = 106 MeV), which is close to the value originally obtained in |27| . 
The order parameter and the susceptibility around this CEP are shown in Fig. |5J 

At /i = 0, as discussed in Introduction, the deconfinement CEP is important as well as the chiral CEP. [In the present 



two-flavor case the chiral CEP trivially lies at m q = 0.] We found the deconfinement CEP at 



788 MeV, T D 



257 MeV) in our model. This value is consistent with the lattice observation in the two-flavor case |l7j and also 
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FIG. 3: The order parameter and susceptibility around the deconfinement critical end-point with mj = 788 MeV and fi = MeV. 



in good agreement with the Gocksch-Ogilvie model 13] . The order parameter and the susceptibility around the 
deconfinement CEP are shown in Fig. |3| Since the deconfinement transition is of second-order, the Polyakov loop 
susceptibility diverges at T = To- In Figs. and|31 we can see that the chiral susceptibility and the Polyakov loop 
susceptibility both have a singularity due to the mixing. 

Discussion We shall briefly summarize characteristic features of our model defined by Eqs. Q and JBJ) below. 

1. The coupling between the Polyakov loop and the chiral condensate is determined uniquely and is consistent with 
the conventional form; in the leading order of the hopping parameter (k) expansion on the lattice, the coupling 
term takes the form of (2k) n -1 ~ Ze" M / T [HEil (see the coupling term of ©). 

2. Because the coupling term is ~ le~ M / 7 ', it goes to zero as T — » 0. Actually the zero temperature system can 
be different from the system at infinitesimally small but finite temperature. At zero temperature, the canonical 
description with the quark triality fixed at zero is likely to be valid 0- 

3. Contrary to naive expectation, the Polyakov loop behavior hardly reflects the singularity associated with the 
chiral phase transition of second-order. [The gross feature of Fig. Q] is hardly changed when m q — except that 
Xg's peak becomes divergent.] This is because the coupling C q i (amplitude between the Polyakov loop and the 
chiral condensate) is proportional to the constituent quark mass, M, and vanishes at the chiral phase transition 
of second-order. Nevertheless, our idea leads to the almost simultaneous crossovers in a robust way. 

4. Figure |21 is an interesting prediction from our model at finite density. Our idea would not necessarily give the 
simultaneous crossovers at high density where Te is too lower than Tq and the Polyakov loop has a long tail. 
This double-peak structure with a sharp peak and a broad bump would be a realistic possibility to be seen in 
the future lattice simulation at high density. 

5. As discussed in the section "Idea", our idea would not work for large m q because of explicit symmetry breaking. 
In general the pseudo-critical temperature T x gets larger with increasing m q . Actually the chiral susceptibility 
peak calculated in the standard NJL model yields T x = 270 MeV(~ T D ) at m q = 167 MeV. For m q > 167 MeV, 
our model leads to a double-peak structure with Ti < T x , that is not prohibited by our dynamical mechanism. 
The chiral susceptibility in Fig. [21 has a second broad bump at much higher temperature than shown in the 
figure. Although the present model embodies our idea and describes the simultaneous crossovers for small m q , 
there must be some further constraints, in particular to impose T^ > T x . Such a condition would realize a 
perfect locking for small m q and cure the failure for larger m q . 

Summary We proposed an idea to realize the simultaneous crossovers of deconfinement and chiral restoration, 
which turned out to work well for small quark mass. We demonstrated the idea by using a chiral effective model with 
the Polyakov loop. The model study yields the chiral CEP at (/ie = 321 MeV, Te = 106 MeV) and the deconfinement 
CEP at (m q = 788 MeV, Tp = 257 MeV). The Polyakov loop susceptibility and the chiral susceptibility both diverge 
at the CEP due to the mixing effect. 

Since the thermodynamic potential in our model is a function of M 2 (M being the constituent quark mass), the 
mixing effect oc M is small at T = Tq. Then our idea plays an essential role to attract one crossover to the other. 
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Also we presented a prediction from our model at finite baryon density. At sufficiently high density we can expect 
that the Polyakov loop has a double-peak structure with a sharp peak from the mixing and a broad bump from a 
remnant of deconfinement. 

The present model lacks some mechanism necessary to sustain the locking for m q > 167 MeV. In other words, this 
result suggests that some dynamical mechanism is further needed in order that the chiral and deconfinement CEPs 
are connected by a single crossover line |Tfjj . The scenario of [lfij requires something beyond the mixing argument 
and the present idea to lock two phenomena. Although this is still an open question, we believe that our model can 
contain correct physics at least for small m q and can be a simple starting point to examine the underlying relation 
between deconfinement and chiral restoration not only at finite temperature but at finite baryon density also. 
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